Numerical Analysis and Scientific Computing Preprint Seria An adaptive octree finite element method for PDEs posed on surfaces

نویسندگان

  • A. Chernyshenko
  • M. A. Olshanskii
  • Alexey Y. Chernyshenko
  • Maxim A. Olshanskii
چکیده

The paper develops a finite element method for partial differential equations posed on hypersurfaces in R , N = 2, 3. The method uses traces of bulk finite element functions on a surface embedded in a volumetric domain. The bulk finite element space is defined on an octree grid which is locally refined or coarsened depending on error indicators and estimated values of the surface curvatures. The cartesian structure of the bulk mesh leads to easy and efficient adaptation process, while the trace finite element method makes fitting the mesh to the surface unnecessary. The number of degrees of freedom involved in computations is consistent with the two-dimension nature of surface PDEs. No parametrization of the surface is required; it can be given implicitly by a level set function. In practice, a variant of the marching cubes method is used to recover the surface with the second order accuracy. We prove the optimal order of accuracy for the trace finite element method in H and L surface norms for a problem with smooth solution and quasi-uniform mesh refinement. Experiments with less regular problems demonstrate optimal convergence with respect to the number of degrees of freedom, if grid adaptation is based on an appropriate error indicator. The paper shows results of numerical experiments for a variety of geometries and problems, including advection-diffusion equations on surfaces. Analysis and numerical results of the paper suggest that combination of cartesian adaptive meshes and the unfitted (trace) finite elements provide simple, efficient, and reliable tool for numerical treatment of PDEs posed on surfaces.

برای دانلود رایگان متن کامل این مقاله و بیش از 32 میلیون مقاله دیگر ابتدا ثبت نام کنید

ثبت نام

اگر عضو سایت هستید لطفا وارد حساب کاربری خود شوید

منابع مشابه

Numerical Analysis and Scientific Computing Preprint Seria Error analysis of a space-time finite element method for solving PDEs on evolving surfaces

In this paper we present an error analysis of an Eulerian finite element method for solving parabolic partial differential equations posed on evolving hypersurfaces in Rd, d = 2, 3. The method employs discontinuous piecewise linear in time – continuous piecewise linear in space finite elements and is based on a space-time weak formulation of a surface PDE problem. Trial and test surface finite ...

متن کامل

Numerical Analysis and Scientific Computing Preprint Seria A narrow-band unfitted finite element method for elliptic PDEs posed on surfaces

The paper studies a method for solving elliptic partial differential equations posed on hypersurfaces in RN , N = 2, 3. The method allows a surface to be given implicitly as a zero level of a level set function. A surface equation is extended to a narrow-band neighborhood of the surface. The resulting extended equation is a non-degenerate PDE and it is solved on a bulk mesh that is unaligned to...

متن کامل

Numerical Analysis and Scientific Computing Preprint Seria Non-degenerate Eulerian finite element method for solving PDEs on surfaces

The paper studies a method for solving elliptic partial differential equations posed on hypersurfaces in RN , N = 2, 3. The method builds upon the formulation introduced in Bertalmio et al., J. Comput. Phys., 174 (2001), 759–780., where a surface equation is extended to a neighborhood of the surface. The resulting degenerate PDE is then solved in one dimension higher, but can be solved on a mes...

متن کامل

Numerical Analysis and Scientific Computing Preprint Seria Numerical integration over implicitly defined domains for higher order unfitted finite element methods

The paper studies several approaches to numerical integration over a domain defined implicitly by an indicator function such as the level set function. The integration methods are based on subdivision, moment–fitting, local quasi-parametrization and Monte-Carlo techniques. As an application of these techniques, the paper addresses numerical solution of elliptic PDEs posed on volumetric domains ...

متن کامل

ذخیره در منابع من


  با ذخیره ی این منبع در منابع من، دسترسی به آن را برای استفاده های بعدی آسان تر کنید

عنوان ژورنال:

دوره   شماره 

صفحات  -

تاریخ انتشار 2014